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ABSTRACT 


This paper is concerned with optimal control problems for the von Karman equations with 
distributed controls. We first show that optimal solutions exist. We then show that Lagrange 
multipliers may be used to enforce the constraints and derive an optimality system from which 
optimal states and controls may be deduced. Finally we define finite element approximations of 
solutions for the optimality system and derive error estimates for the approximations. 
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1. Introduction 


The von Karman equations for a clamped plate are given by (see, e.g., [7] or [9]) 


and 

where 


A 2 V»i + = 0 in ft 


A 2 V >2 - [V>i > i’i] = A g in ft 
’ dx\ dx\ dx\ dxj dx\X 2 dx\X 2 


Here, ft is a bounded, convex polygonal domain in 1R 2 , V >1 denotes the Airy stress function, ip 2 
denotes the deflection of the plate in the direction normal to the plate, and A g is an external load 
normal to the plate which depends on the loading parameter A. 

The boundary conditions on T = dfl are given as follow: 

, din , di>2 ft ^ 

* = *r=* = *r = 0 °" r ’ 

where d(*)/dn denotes the normal derivative in the direction of the outer normal to I\ 

For reasons to be explained later we introduce appropriate rescalings, i.e., by replacing tp i by 
Xipi, ip 2 by X / tp2 > and g by Xg. Then we can rewrite the Karman equations as follows: 



A 2 V»i + ^[^ 2 , in] = 0 in ft , 

(1.1) 

and 

A 2 $j - A [-01 , -02] = A gr in ft, 

(1.2) 

, din di>2 n 7-1 

= dn = *= dn = ° °" r - 

(1.3) 


We introduce the standard Sobolev spaces (see, e.g., [1]) 


(ft) = { 


i> g 2 (ft) 


^ = 0,^ = 0 on r 
on 


}• 


ff 0 2 (ft)= (tf 0 2 (ft)) 2 , ff- 2 (ft)=(ff 2 (ft))*, and ff~ 2 (ft) = (ff 0 2 (ft))* , 

the bilinear form 

a{i>,<j>) = I AipAcfxlSl ff 2 (ft) 

J n 


and the L 2 (ft)-inner product 


(if,4>)= t i>i>dSl V L 2 (Cl ) . 
J o 


Then we may define the following weak formulation of the von Karman equations (1.1)-(1.3): find 
i> = (i}\,xp 2 ) T G H 2 (Sl) such that 

a (i>i,4>i) + ^([^2,^2],^!) = 0 V <£1 G Hq(SI) (1.4) 
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and 


= \(g,fo) VfoeHSW). (1.5) 

Here the notation (•, •) stands for the L 2 (fl)-inner product and (•, •) the duality pairing. Using the 
identity 

([tM,C) = (ty.CU) (1.6) 

one can show that for each g 6 (1.4)-(1.5) possesses at least one solution i/> = (V’i ? V , 2) T € 

Hq(SI) and that all solutions of (1.4)-(1.5) satisfy the a priori estimate 

HV'i lh 4* ||^2||2 < C'||p||_2 ; (1.7) 

see, e.g., [9], for details. In the sequel a solution to (1.1)-(1.3) will be understood in the sense of 
(1.4)-(1.5). T 

Given a desired state ij ) 0 = (V'io^o) € L 2 (ft), we define the functional 

+ ~ ^ 20 ^ 2 ) ^ j g 2 dSl ^ ^ 

for all V* = (V’IjV^) 7 ' G tf 0 2 (ft) and g G L 2 (fi). We then consider the following optimal control 
problem associated with the von Karman plate equations: 

min { J(^,g) | 5 J6 0} subject to (1.4)-(1.5), (1.9) 

where 0 is a convex, closed subset of L 2 (f2). The physical interpretation of this optimal control 
problem is that we wish to match a desired Airy stress distribution and a deflection distribution 
on the entire plate by choosing an (optimal) external load g from an admissible set 0. 

Our plan is as follows. In §2, we show that an optimal solution exists. In §3, we use Lagrange 
multiplier techniques to derive an optimality system from which optimal states and controls can be 
deduced. We then specialize to the case 0 = L 2 (fi) and derive an optimality system of equations. 
In §4, we define conforming finite element approximations of the optimality system of equations 
and derive error estimates. 


2. Existence of an optimal solution 

We will prove in this section that an optimal solution exsits for the minimization problem 

(1.9). 

Theorem 2.1. There exists a ( <f>,g ) 6 x 0 such that (1.8) is minimized subject to 

(1.4)-(1.5). 

Proof. We choose a £f(°) G 0 and let ij)^ = € Hq(CI) be a solution of 

and 

a(4°\<f>2)- A([W 0) ,4 0) ]^2) = Ate<°>,&) HfcSl). 
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The existence of such a V» (0) was established in [9]. We see that (^ (0) ,S (0) ) € ^(fl)x0 satisfies the 
constraints (1.4)-( 1.5) and J(V» (0) ,5 (0) ) < oo. We also note that g) > 0 for all (i/>, g) £ X x0. 
Thus we may choose a minimizing sequence {(V» (n) ,fif( n ))} C tf^(fl) x 0 such that for some constant 
M > 0, 

J{j> (n) ,9 (n) ) < M, (2.1) 

<4 n \fo) + ^{[ti n \4 n) Ui) = 0 V*e JT 0 2 (fl) (2.2) 

and 

a(4 n) ,fo) ~ M[W n) ,V4 n) ],<k) = \(gW,<h) V <f> 2 € JTftfl). (2.3) 

Using (2.1) and the definition of the functional J - we deduce that the sequence {||srt n > ||o} is bounded 
so that by the a priori estimate (1.7), the sequence is bounded as well. Hence by choosing 

a subsequence, we have that 

g(n) g j n 

and 

^ (n) — in tf 2 (11) 

for some g £ L 2 (H) and i[> £ Hq(SI). Using the convexity and closedness of 0, we deduce that 
g £ 0. By the compact imbedding tf 2 ( 11) <-»<-»• L 2 (ll), we obtain V» (n) -> i/> in L 2 (H) so that for 
each <f> = (<f>i,(f> 2 ) T £ Hq(Q,), we have that (also using the continuous imbedding Hq(Q,) t — > L°°( 11)) 

n^Joo \ ([^ 2 Tl) ’^ 2 n) ])^i) = Jirn i ([ip ( 2 n \<t>i\, 4 n) ) 

= ^([^2,<Ai]> fh) = ^([V»2,^2],^l) 

and 

Jim - ([V4 n) ,^ n) ],^) = Km - ([^* ) ,^a],^ B) ) 

= ~([V’l»^2]»V > 2) = -([V’i,^],^) • 

Hence we may pass to the limit in (2.2)-(2.3) to show that (lf>,g) satisfies (1.4)-(1.5). 

Finally, we use the sequential lower semi- continuity of the functional J to obtain that 

J(S>,g) < Urn J(V» (n) , 5 (n) ) . 

n _ ►oo v ' 

Hence, we have proved that a solution (rj>, g) £ tf 0 2 (l 1) x 0 exists that minimizes (1.8) subject to 
(1.4)-(1.5). □ 


3. An optimality system of equations 

In this section we assume that x 0 is an optimal solution for the minimization 

problem (1.9) and we attemp to characterize the optimal solution as the solution for a system 
of partial differential equations. To be precise, we use Lagrange multiplier rules to derive an 
optimality system of equations. 
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We define the Lagrangian for the constrained minimization problem (1.9) as follows: 


+ a(V>2,»?2) - A([- 0 i, ^2 ], t? 2 ) - Kg,m)} 

for all W,g,ri) 6 H$(£l) x L 2 (ft) x 

By formally taking variations in the Lagrangian with respect to \j) and g , we obtain: 

a(Ci,»n)-M[^,i* 2 ],Ci) = A (^1 -^io,Ci) VCi 6 ff 0 3 («). (3-1) 

a(C2,»?2) + A([^2, 771], C2) - A([V’1, 772], C2) 2 

= A(V>2-V’20,C 2 ) V C2 € Hq(Q) ' 

and A A 

-(z,z) + X(z,rj2)- - (g,g)~ A(5,t? 2 ) > 0 Vz€ 0 . ( 3 . 3 ) 

For each € € ( 0 , 1 ) and each i G 0 , we set 2 = ef + (1 — e)g € 0 in ( 3 . 3 ) to obtain 

y (t- g, t-g) + 6 (t- g,g)+ e(t- g, 772) >0 Vt 6 0 
so that, after dividing by c > 0 and then letting c — ► 0+, we obtain 

(t-g,g + v 2 ) > 0 v t € 0 . (3.4) 


Now we show that there does exist an 17 € Hq(SI) satisfying ( 3 . 1 )-( 3 . 3 ), or equivalently, ( 3 . 1 )- 
( 3 . 2 ) and ( 3 . 4 ), so that we are justified to compute a triplet (V>,<7,»7) from ( 1 . 4 )-( 1 . 5 ), ( 3.1 )-( 3 . 2 ) 
and ( 3 . 4 ). In this paper we will not address the uniqueness of solutions for the system ( 1 . 4 )-( 1 . 5 ), 
( 3 . 1 )-( 3 . 2 ) and ( 3 . 4 ). 

We first quote the following abstract Lagrange multiplier rule whose proof can be found in 

[1 °]. 

Theorem 3.1. Let Xi and Xi be two Banach spaces and 0 an arbitrary set. Suppose J is a 
functional on X 1 x 0 and M a mapping from X\ x 0 to Xi. Assume that ( u,g ) is a solution to 
the following constrained minimization problem: 

M(u,g ) = 0 and there exists an e > 0 such that J(u,g) < J(v,z) 

( 3 . 5 ) 

for all ( v,z ) such that ||u — u||^j < e and M(y,z) = 0 . 

Let U be an open neighborhood of u in X\ . Assume further that the following conditions are 
satisfied: 

(A) for each 2 6 0 , vh ff(y,z) and v >-+ M(v,z ) are Frechet-differentiable at v = u\ 

(B) for any v G U, 21,22 € 0 , and 7 6 [ 0 , 1 ], there exists a 2 7 = z^(v,zi,Z2) such that 

A/(v, 2 7 ) = 7M(u,2l) + (1 - 7)Af(u,22) 
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and 

j{v,z^) < lJ{v,zi) + {\-i)J(v,z 2 )\ 

(C) Range(M u (u, g)) is closed with a finite codimension; 
and 

(D) the algeraic sum M u (u,g)X\ + M(u, 0) contains 0 € X 2 as an interior point, 

where M u (u,g) denotes the Frechet derivative of M with respect to u. Then, there exists a a rj € X% 
such that 

(r],Mu(u,g)v)- (J u (u,g),v) = 0 V v € X\ 

(or equivelently, [ M u (u,g)]*r] - J u (u,g ) = 0 ) and 

minC(u,z,rj) = C(u,g,rj), 

where £(u, < 7 , 7 /) = J(u,g) — (r i /, M(u,g)} is the Lagrangian for the constrained minimization problem 

(2.5) and where Ju(u^g) denotes the Frechet derivative of J with respect to u. □ 

We now define some spaces and operators in order to rewrite the constraint equations (1.4)- 

(1.5) into a form that will facilitate our verification of the assumptions in Theorem 3.1. 

We define the spaces X = Y = /7 _ 2 (fi), G = L 2 (ft), and Z = L l (tl). By compact 

imbedding results, Z Y . The control set 0 is a closed, convex subset of G = L 2 (fl). We also 
note that Y = X *. 

We define the continuous linear operator A £ C{X]Y) as follows: for if) G X = /?o(fl), 
Ail) = f £ Y = H ~ 2 (fi) if and only if 

a(1n,4n) = (fu<h) V fa e H$(Sl) 

and 

a(tl>2,fo) = (f2,fo) 

It can be easily verified that A is self-adjoint, invertible, and A~ l 6 C(Y;X). 

We define the (differentiable) nonlinear mapping N : X — ► Y by 

w) =(-!*;«) 


or equivalently 


(N(4>),4>) = 1 ([^ 2 ,^ 2 ], * 1 ) - ([V>i , fc], fc) V * = (<£, , € X 


and define K : g € L 2 (fl) 


or equivalently, 


**-(!)• 

(Kg,#) = -(g,<f>2) v^ = (^)gx. 
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Clearly, the constraint equations (1.4)-(1.5) can be expressed as 

Arf> + A N(rj>) + A Kg = 0 . 

We are now in a position to prove the existence of a Lagrange multiplier t) 6 H$ (fl) that 
satisfies (3.1)-(3.2) and (3.4). 

Theorem 3.2. Assume G Hq(SI)xQ is an optimal solution that minimizes (1.8) subject 
to (1.4)-(1.5). Then, for almost all A € A, there exists a Lagrange multiplier ij £ Hq(£ 1) satisfying 
the Euler equations (3.1 )-(3.2) and (3.4). 

Proof. We define a mapping M : x @ — > ff -2 (ft) (i.e., X x 0 — ► Y) by 

M{4>,z) = A<f> + \N(<l>) + \Kz V {<f>,z) £ x 0. 

We see that minimizing (1.8) subject to (1.4)-(1.5) can be stated as minimizing J(iJ),g) subject to 
M(if>,g) = 0. We first verify that the hypotheses (A)-(C) of Theorem 3.1 hold with X\ = X and 
X 2 = Y. 

Obviously, for each z G 0, the mappings <f> ■-> J(<f>,z) and <f> <->■ M(<f>,z) are both Frechet 
differentiable, i.e., (A) holds. 

Since 0 is convex and the mapping K is linear, we have that for any <f> £ Hq(£1), z\ £ 0, 
z 2 6 0 and 7 £ (0, 1), the element z = jzi + (1 — 7)^2 belongs to 0 and 

Zy) = A<f> + A + A Kz-t 

= 7 (A<f> + A N(4>)) + (1 - 7) {A<t> + A 
+ \{yKzi + (1 - 7)^2)) 

= 7M(^,zi) + (1 - 7)Af(^,z 2 ) • 

Moreover, the mapping g ►-+ j | £7 1 1 5 is obviously convex from 0 into IR so that 
J(<f>,z-,) = ^ ||^- ^ollo + ^ lll^i + (1 - 7)^2||§ 

< 5(7 U - 4,11$ + (i - 7 ) 114- 4,113) + \ (7 INIS + (i - 7 ) INIS) 

= + - 7 ) J(<t>,Z 2 )- 

Thus, (B) holds. 

The operator Af^,(^, g) from X to Y is defined by 

AfyW, 9) • 4> = A<t> + A W{4>) -<j> V<f>£X = H>(Sl), 

or simply, 

M^,g)= A + \N'W), 
where for any ^ 6 X, the operator : X -* Y is given by 

(-[*.« **,*]) 
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Thus, using the definition of [*,•], we obtain that iW(^>) • <f> 6 = Z. Using the compact 

imbedding Z <-><-*■ V , we see that N'(iJ)) is a compact operator from I to I. As a result, 
= A + A = A(l + \ A - * N'(ip)) is a Fredholm operator so that it has a closed 

range with a finite codimension, i.e., (C) holds. 

We now verify that the hypothesis (D) holds for almost all A. In fact, if A ^ 0 and (1/A) £ 
^(-A-'N'W)), where (r(-.4 _1 N'(i/))) denotes the spectrum of the operator (-A -1 N'(ip)), then 
it follows that 

Y = Range(A(I + A A -1 N'(tf>))) = Range(M r p(iJ),g)) 

so that Range(M 1 j ) (il),g )) contains 0 € Y as an interior point, i.e., (D) holds. Since the spectrum 
of a compact operator is at most countable, we conclude that for almost all A, (D) holds. 

Hence, by Theorem 3.1, we obtain that for almost all A, there exists a T) £ Y* = Hq(£1) such 

that 

A*tl+\[N'{1>)]*-ii-JfW,g) = 0 (3.6) 

and 

^,9,v) < £(^,z,» 7) Vz£&. (3.7) 

Recalling the definition of the operators A and N'{xj)) as well as the fact that A is self-adjoint, i.e., 
A* = A, we see that (3.6) is equivalent to (3.1)-(3.2). Using the definition of the Lagrangian, we 
easily see that (3.7) is simply a restatement of (3.3), the latter being equivalent to (3.4). Thus we 
have proved that a Lagrange multiplier tj 6 ffg(il) exists such that (3.1)-(3.2) and (3.4) hold. □ 

Remark. Now we can explain the rescalings that were introduced in §1. With these rescalings, 
the parameter A appears in the form of an eigen value in the operator M^{if),g) = A + A 
so that spectrum theories for compact operators can be readily employed in the proof of Theorem 
3.2. □ 

So far 0 has only been assumed to be a closed, convex subset of L 2 (Sl). In the sequel we 
specialize to the case 0 = L 2 (?l). Then we can obtain two improvements in Theorem 3.2, the first 
being that the results of Theorem 3.2 now hold for all A, the second being that (3.4) becomes an 
equality. 

Theorem 3.3. Assume 0 = L 2 (fl) and (i>,g) € Hq(£1) x L 2 (Q,) is an optimal solution that 
minimizes (1.8) subject to (1.4)-(1.5). Then, for all A £ A, there exists a Lagrange multiplier 
f) £ Hq(SI) satisfying the Euler equations (3.1)-(3.2) and 

(t,9 + i 72) = 0 VteT 2 (ft). (3.8) 

Proof. In view of Theorem 3.2 and its proof, we only need to verify (D) when (1/A) e 
A -1 N'(u )) , i.e., we need to show that for each / £ H~ 2 (tl), there exists a g £ L 2 (Sl) and a 
£ H$(Cl) such that 

Ai> + \N'(iP).j> + \Kg = f, 

or equivalently, 

+ H 2 (Q) (3.9) 
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and 


(3.10) 


a(i>2,<h) - A([$i,tfe], fa) - - Kg,h) 

= (h,<h) v ^>2 g Him) . 

To show this, we first let i[> 6 Him) be a solution of 

a(i>i,4>i) + x{[tf> 2 ,i> 2 ],<l>i) = (/i , 4>\ ) v <^i e Him) 

and 

a(i>2,(h) ~ = (f2,<h) V 4> 2 g Him). 

The existence of such a if) can be shown in a manner similar to that for showing the existence of 
a solution to the von Karman equation; the key step is that by adding the two equations with the 
test function <f> replaced by if), we have the a priori estimate 

a(^>$i) + a (^2,&) = </i,V>i) + </2,V>2). 

Now, having chosen such a rf) y we simply set g = —[^ 1 ,^ 2 ]. Note that regularity results for the 
biharmonic equation applied to (1.4)-(1.5) yield if) £ 77 4 (fi) (see [3]). Hence, using imbedding 
theorems we deduce that g £ L 2 (Q). It is obvious that g and if) satisfy (3.9)-(3.10). Thus, we have 
verified (D) so that by Theorem 3.1, there exists a r) £ X such that (3.1)-(3.2) and (3.4) hold. 
Finally, (3.4) trivially reduces to (3.8) in the present case of 0 = L 2 (ft). □ 

Combining the results of §2 and this section we see that we have proved the existence of a 
triplet (if>,g,r)) £ H$(£l) x L 2 (fi) X H$(Q,) suc ^ that the system (1.4)-(1.5), (3.1)-(3.2) and (3.8) 
is satisfied. For convenience we collect these equations here to form an Optimality System of 


Equations: 

o(ih,4n) + ^([l>2M,<h) =0 v </>i g Him ) , (3.11) 

a(^ 2 ,</> 2 ) - A([^i,V’ 2 ],<^ 2 ) = >i(g,fa) V(^Gffo(fl), (3.12) 

a(Cum) - A«*,w],Ci) = A(tfi - tfio, Cl) VCi G Him), (3-13) 

fl(C2,7?2) + A([t^2, r?l], C2) — A([V>1, »72], C2) 

= A(V>2-V>20,C2) V C2 G Him) { ' 

and 

(*,<7 + f/2) = 0 VtGL 2 (fi). (3.15) 


4. Finite element approximations and error estimates 

4.1. Definition of finite element approximations, 

A finite element discretization of the optimality system (3.11)-(3.15) is defined in the usual 
manner. For simplicity, we will only study conforming finite element approximations in this paper. 
However, the error estimation techniques used in this paper are equally applicable to mixed finite 
element approximations based on the Hellan-Hermann-Johnson scheme for biharmonic equations. 
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See [4] (also [5]) for the definition and discussions of the mixed Hellan-Hermann-Johnson scheme 
for biharmonic equations. 

We first choose families of finite dimensional subspaces X h C Hq(£1) and G h C Z» 2 (ft) param- 
eterized by a parameter h that tends to zero and satisfying the following approximation properties: 
there exists a constant C and an integer r such that 

inf U -<!>%< Ch^\\(f>\\ m+2 , V <f> £ H m+2 (Sl ) , 1 < m < r (4.1) 

<t> h ex h 

and 

inf \\z - z h \\o < Ch m \\z\\ m , V z £ tf m (ft) , 1 < m < r . (4.2) 

z h £G h 

One may consult, e.g., [2] and [6] for some finite element spaces satisfying (4.1) and (4.2). For 
example, one may choose X h = V h x V h where V h is the piecewise quintic-C 1 (ft) finite element 
space constrained to satisfy the given boundary conditions and defined with respect to a family of 
triangulations of ft. In this case, h is a measure of the grid size. 

Once the approximating spaces have been chosen, we may formulate the approximate problem 
for the optimality system (3.1 1)-(3.15): seek tj) h £ X h , g h £ G h , and f) h £ X h such that 


<t>\) + \ rt], <t>i) = 0 V#€K\ (4.3) 

= Ka h ,4>2 ) v $ £ v h , (4.4) 

«(C 1, Vi) - A ([V »2 > V2I Cl) = Ki’i - -010 , Cl) V Cl 6 v h , (4.5) 

«(C 2 \ V 2 ) + HM, 17 ?], C^) - A(K, % fc ], C 2 ") = Ki >2 - i> 20 ,cl l ) V e v h (4.6) 

and 

(z k ,g h + 773 ) = 0 V z h £ G h . (4.7) 


Note that in the last equation, if G h = V h , then we have g h = —rj^ so that the variable g h can be 
eliminated to simplify the approximate problem. But in general we have to deal with the entire 
system (4.3)-(4.7). 

4.2. Quotation of results concerning the approximation of a class of nonlinear prob- 
lems. 

The error estimate to be derived in Section 4.3 makes use of results of [4] and [8] concerning the 
approximation of a class of nonlinear problems. These results imply that, under certain hypotheses, 
the error of approximation of solutions of certain nonlinear problems is basically the same as the 
error of approximation of solutions of related linear problems. Here, for the sake of completeness, 
we will state the relevant results, specialized to our needs. 

The nonlinear problems considered in [4] and [8] are of the following type. For given A £ A, 
we seek ip £ X such that 

ft(A,^) = ^ + r<7(A,VO = 0, (4.8) 

where T £ £(y;X), Q is a C 2 mapping from A x X into X and y are Banach spaces, and A 
is a compact interval of 1R. We say that {(A,^(A)) : A £ A} is a branch of solutions of (4.8) if 
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A — *• V’(A) is a continuous function from A into X such that H( A, tf>( A)) = 0. The branch is called 
a regular branch if we also have that 7f^(A, V’(A)) is an isomorphism from X into X for all A € A. 
Here, %,/,(•,■) denotes the Frechet derivative of ?f(*, •) with respect to the second argument. We 
assume that there exists another Banach space Z , contained in with continuous imbedding, 
such that 

MA ,4))e£(X;Z) V A € A and ip € X , (4.9) 

where £/,/,(•,•) denotes the Frechet derivative of £/(•,•) with respect to the second argument. 

Approximations are defined by introducing a subspace X h C X and an approximating operator 
T h e C(y;X h ). Then, given A 6 A, we seek ip h 6 X h such that 

H h (\, rp h ) = il) h + T h Q(\, ip h ) = 0 . (4.10) 

Concerning the operator T h , we assume the approximation properties 

lim \\(T h - T)u\\ x = 0 V u e y (4.11) 

h — ►O 

and 

«>n||(r*-r)|| r( *,* ) = o. (4.i2) 

Note that whenever the imbedding Z C y is compact, (4.12) follows from (4.11) and, moreover, 

(4.9) implies that the operator TGip(\ y ip) E C(X\X) is compact. 

We can now state the result of [4] or [8] that will be used in the sequel. In the statement of 
the theorem, D 2 Q represents any and all second Frechet derivatives of Q. 

Theorem 4.1. Let X and y be Banach spaces and A a compact subset of JR,. Assume that 
Q is a C 2 mapping from A X X into y and that D 2 Q is bounded on all bounded sets of A X X. 
Assume that (4.9), (4.11), and (4.12) hold and that {i/>(A);A E A} is a branch of regular solutions 
of ( 4.8). Then, there exists a neighborhood O of the origin in X and , for h < ho small enough , a 
unique C 2 function A — » ^(A) E X h such that {^(A); A E A} is a branch of regular solutions of 

(4.10) and / 0 /l ( A) — 'ip(X) E O for all A E A. Moreover , there exists a constant C > 0, independent 
of h and A, such that 

||^(A)-^(A)||Ar<C||(r*-r)f?(A,^(A))|U VA € A. □ (4.13) 


4.3. Error estimates. 

In this section we will derive error estimates for the finite element approximations of solutions 
of the optimality system. 

We begin by fitting the optimality system and its finite element approximations into the 
abstract framework in §4.2. 

We set X = X x G x Y*, y = Y x X*, Z = Z x Z, and X h = X h x G h x (Y*) h . Then 
we have that Z is compactly imbedded into y. (We recall that the spaces X = Hq(Q) = Y*, 
Y = H~ 2 (£l) = X*, Z = IA(fl) and G = T 2 (fl) were all introduced in §3.) We define T € C{y\X) 
as follows: T(f, f) = (V>,5,^) for (r,f) € y and € X if and only if 

«(&*) = (*,*) (4.14) 
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and 


a(i),w) = (f,u>) Vw6 H$(Sl) 


(4.15) 


(t,g + i 72) = 0 Vt€L 2 (Q). (4.16) 

Similarly, the operator T h 6 C(y;X h ) is defined as follows: T h (r,f) = ,g h ,f) h ) for (r,f) € y 
and ( rf> ,g h ,ij h ) e X h if and only if 

a& h A h ) = (r,4> h ) (4.17) 

a(fj h ,u h ) = (f,u h ) Vu h eX h (4.18) 

and 

(*\$ fc + $) = 0 V t h £. G h . (4.19) 

The nonlinear mapping £7 : A x X -» is defined as follows: £7 (A ,(r/),g,ij)) = (r,r) for A £ A, 
{$■> 91 * 1 ) £ X, and (r,f) € 3^ if and only if 

f=( ) 

U*,*J + A gj 

and 

T=( - A ^ 2 ’ ^ t A (^ ” ^i°) ^ 

\ — A[^2 ,^i] + A[V> 1 , 172] -f A(l/>2 — ^20)/ ’ 

or equivalently, 

(a(A,(^5^)),(^,c)) 

= _ 2 ([^2,^2],<Ai) + A ([^1, i> 2 ],<h) + A (g,<t>2) 

+ A([^2,772],Cl) + A (V*1 - V’NbCl) - A ([V>2, ffl ],&) 

+ A([&,HCa) + A(& - ifeo.Ca) V (<f>,Q € ^o 2 («) x ff 0 2 (fl)- 

Using the operators N and K introduced in §3, we can also write the defining equations for 
(7 (A ,(^,g,ij)) as follows: 

f = \N$) + \Kg (4.20) 

and 

f = A [N'(ij))]*r) - A (^ - ^0) • (4.21) 

It is easy to see that the optimality system (3.11)-(3.15) and its finite dimensional counterpart 
(4.3)-(4.7) can be written as 

(&$»*?) + r0(A,(lfc £,»?)) =0 

and 

(^,^,t? /l ) + W(A,(/,^,^)) = 0, 

respectively, i.e., with V> = (V»,5>»?) and ^ = (J> h ,g h ,i) h ), in the form of (4.8) and (4.10), respec- 
tively. 
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We now turn to the verification of the assumptions in Theorem 4.1. We first examine the 
approximation results for the linear operators T and T h . 

Lemma 4.2. For every (f,f) G y, 

/<m||(T-r»)(f,f)|U = 0. 

Furthermore, ifT(f,f) = (^,</,»j) G /f m+2 (ft) x x H m+2 ($l), then 

II (T - T h )(r,f)\\x < Chm (|W| ro+2 + ||sf||m 4- II^IU+ 2 ) . 

Proof. For each given (f, f) 6 y , let {^>,g,fj) = T(f,f) and (V> ,g h ,t) = T h {r,f), i.e., 
(Vb5>*?) € /V and {if ,g h ,t) G X h satisfy (4.14)-(4.16) and (4.17)-(4.19), respectively. From the 
well-known results concerning the approximation of biharmonic equations (see, e.g., [2] and [6]), 
we obtain that 

lim ||V> — = 0 

h-*o" " 

and 

lim\\fi-f,%=0. (4.22) 

A— ► 0 

Furthermore, if ip G // m+2 (n) and i/ € /f m + 2 (fl), then 

||^-^||2 < Ch m \\xj)\\ m+2 

and 

||^-^||2 < C7l m ||»?||m+2 . 

(4.16) and (4.19) yields ( g — g h ,z h ) = —(772 — V 2 i zh ) f° r all z h G G h . We then have that 
ll<7-<7' l llo = {9-9 h ,9~9 h ) 

= (9 ~ <j h ,9 ~ z h ) + (§ ~ 9 h ,z h - g h ) 

= (g -g h ,g- z h ) - {in - V2’ zh ~9 h ) 

<lls-^IIS + jl|j-j‘ll§ + i||w-« 

+ 5 lk k - silo + II* -’) 2 l |lo + jlls-S'llo 

so that 

II 9 - S fc ||o < C (||<7 - Z*||o + ||*/2 - V 2 ||o) . 

It follows from (4.2) and (4.22) that lim ^ 0 ||<7 — fif^llo = 0. Hence, 

lim || (T - T*)(f, f )||* = Um (||$ - * k || a + ||s - s'llo + ||$ - th) = 0 . 

Furthermore, if (V>,<7,*?) G H m + 2 (£l) x x H m + 2 (Cl), then 

11$ - ^llo < Ch m (\\g\\ m + pllm+ 2 ) 
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so that 

H(r-r‘)(f,f)iu = ||J- *‘|| 2 + ns - ff^Ho + ||« - «*||, 

< cfc”.(iwu +! + ||ffiu + IHHU+O • o 

Next, we examine the derivatives of the mapping Q. 

Lemma 4.3. The mapping Q is twice continuously Frechet differentiable ; the second order 
derivative is bounded on all bounded subsets of Ax X ; and the first order derivative maps X into 
Z continuously, i.e., g €C(X-,Z). 

Proof. Using imbedding theorems, Holder inequalities and the fact that Q consists of poly- 
nomial maps, we deduce that the mapping Q is twice continuously Frechet differentiable and the 
second order derivative is bounded on all bounded subsets of A x X. 

A simple calculation shows that g ^(A, (ip,g,r))) e C(X\y) is given by 




N'W)‘i>+ Kg 
[N"(rf)) ■ ij)]* • T) + [N'(tf))Y 


V-'l> 


») 


where the operators N and K along with their derivatives were introduced in §3. Sobolev imbedding 
theorems imply that 

N ' w -*= (- ( *.fr3L*]) 

g tW ) e*m 

stnd 

«»»»)]••*)•?= « **<«>• 

Of course, ij) G From the definition of the operator K we see that K maps L 2 (fl) into 

{0} x L 2 (Sl) C Hence the operator ^ (A, (V>,<7,77)) maps X into Z. Moreover, using 

Holder inequalities and Sobolev imbedding theorems it is easy to see that 


«•’>)) •(iM>#)llt , <o) 

< CWlfe + Hfflls + 1) (||*||2 + llffllo + llfflb) 


so that 

»(*•(*».»)) €£(*;.Z). 0 


Next, we recall the notion of regular solutions. A solution (ip(\),g(\),r)(\)) of the optimality 
system (3.1 1)-(3.15) is called regular if 


(/ + Tg(il>,g,ri)(\,(il>,g,T)))') -{ip,g,n) = (IAC) 

has a unique solution {^),g,ff) € X for any (|,h,<) € X. As T is invertible, (l,h,Q = T(f,f) with 
(f,f) G y = x H~ 2 { fl). Thus, (^>(A),p(A),t^(A)) is a regular solution of the optimality 

system (3. 1 1)-(3.15) if 


(7 + g, »;)(A, (ip, g,r)))^j • {i>,~g,rj) = T(f, f) 
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has a unique solution (V>,<7,ij) € X for any (f , f) 6 3>, or equivalently, the following system (for 
the unknowns (i>,g,fi)) is uniquely solvable for any (f,f) 6 y-. 

a(v>i,<^i)+ ((^2 ,v>2],<£i) = (n,<f>i) V<£i € #o(ft), 

a(i> 2 ,<f> 2 ) - A([V» ui> 2 ],<h) - A([^i,V»2],^2) - Hg,<t> 2) 

= (n,<M y <f >2 € Hq\ 

a(Ci) j?i) - A([V>2,^2 ],Ci) - A([V’2, >72], Ci) - A(^i,Ci) 

= <ri,Ci) VCi€^ 0 2 ( 

a(c 2,772) + A([V’2,i?i],C2) + A(['02, f?i]? C2) - A ([V’i 1 ^2], C2) 

A([^1, »72], C2) - A(V>2,C2) = (^2,(2) v C2 € H$(£L) 

and 

(M7 + »?2) = 0 VfeL 2 (ft). 

Note that the linear operator appearing on the left hand side of (4.23)-(4.27) is obtained by 
linearizing the optimality system (3.1 1)-(3.15) about {rj),g,rj). 

Lemma 4.4. For almost all X, solutions (V>(A), g(X), rj(\)) of the optimality system (3.11)- 
(3.15) are regular. 

Proof. The system (4.23)-(4.27) can be rewritten as 

(/ + \TS(tl>,g,ri)yij>,g,fi) = T(f,f), (4.28) 

where the linear operator : X — ► y is defined by 

^,g,v) ■ (^,g,v) = jO ( ^g i1}) {X,(il),g,rj/)) -(^ ,g,fi ) 

= ( N'(il>)-]>+Kg 

It was established in Lemma 4.3 that the mapping Q^ g ^(\,{^},g,ri)) is compact from X into y. 
Now, T e C(y-, X), so that (/ + \TS(if),g,r))) is a compact perturbation of the identity operator 
from X to X. Thus, for almost all A, (4.28), or equivalently (4.23)-(4.27), is uniquely solvable, i.e., 
for almost all A, the solution (V>( X),g(X),t }(\ )) of the optimality system (3.1 1)-(3.15) is regular. 
□ 

We are now in a position to derive error estimates. In the following theorem, we will assume 
that the solution (^>(A), 5(A), 17(A)) of the optimality system (3.11)-(3.15) is regular. Lemma 4.4 
guarantees that this is indeed the case for almost all A. Lemmas 4.2 and 4.3 verified all the 
assumptions in Theorem 4.1. Thus we are led to the following error estimates. 

Theorem 4.5. Assume that A is a compact interval 0/IR+ and that there exists a branch 
{(^(A), <7(A), 17(A)) : A € A} of regular solutions of the optimality system (3.1 1)-(3.15) . Assume 
that the finite element spaces X h and G h satisfy the hypotheses (4.1)-(4.2). Then, there exists a 



(fl), 

ft), 


(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 
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^ > 0 and an ho > 0 such that for h < ho, the discrete optimality system (4.3)-(4.7) has a unique 
branch of regular solutions {{^ h (X),g h (X),r) h (X)) : A € A} satisfying 

{|(^> /l (A) - V *(^)||2 + ||<7 h (^) — ?(A)||o + H^HA) - J?(A)|| 2 } < 6 for all A € A. 

Moreover, 

lim{\\if) h (\) - V»(A)|| 2 + H5KA) - 5(A)|| 0 + \\v h (X) - j?(A)|| 2 } = 0 , 
uniformly in A € A. 

If in addition, the solution of the optimality system satisfies (V’(A), < 7 (A),»j(A)) € H m+2 (Q.) x 
H m (£l) x H m + 2 (Cl) for A e A, then there exists a constant C, independent of h, such that 

WW - *‘(A)|| a + HsW - 9‘(A)||o + IMA) - n*(A )||2 
< Cft-»(||0(A)|U +2 + || S (A)||„ + HMA)|| m «) , 

uniformly in A G A. □ 


5. Conclusions 

We studied an optimal control problem for the von Karman equations in this paper. We first 
gave the mathematical statement of the problem and proved the existence of an optimal solution. 
We then applied the Lagrange multiplier rules to derive an optimality system of equations that the 
optimal state must satisfy (the use of Lagrange multiplier rules was justifed). We finally defined 
finite element approximations for the optimality system and derived optimal error estimates. The 
functional we minimized was a tracking functional (i.e., the tracking of the Airy stress function 
and the normal deflection). The control we used was a distributed control (i.e., the external load). 
However, the methods used in this paper apply equally well to optimal control problems with other 
objectives (e.g., minimizing the stress functional in some areas) and/or other types of controls (e.g, 
boundary controls). 
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